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ABSTRACT 



Within the finite temperature TV = 4 strongly coupled super Yang- 
Mills, we compute the dissipative force on an external quark in the pres- 
ence of evenly distributed heavy quark cloud. This is computed holo- 
graphically by constructing the corresponding gravity dual. We study the 
behaviour of this force as a function of the cloud density. Along the way 
we also analyze the stability of the gravity dual for vector and tensor 
perturbations. 
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Collisions of heavy nuclei are believed to produce QGP, a strongly coupled thermal 
state of matter [H [2] . This motivated many researchers to compute QGP observables 
using gauge/gravity correspondence. Though the gauge theory in question is quite 
different from QCD, this correspondence gives us a systematic prescription to com- 
pute various quantities associated with a strongly coupled thermal system [3j. These 
include viscosity, entropy production, transport coefficients to name a few [H El |6]. 
In this respect, an interesting quantity is the dissipative force experienced by an ex- 
ternal heavy quark moving in a thermal plasma [3 [HI |9]. Within the framework of 
gauge/gravity duality, an external heavy quark is modeled by a fundamental string 
attached to the boundary of an AdS black hole. For A/" = 4 super Yang-Mills, the 
end point of the string carry a fundamental SU{N) charge. The string extends along 
the radial direction of the AdS-Schwarzschild metric. This external quark, with a 
mass proportional to the length of the string, loses its energy as the string trails back 
imparting a drag force. 

Our aim in this paper is to carry out a simple calculation of the drag force for 
thermal A/" = 4 super Yang-Mills on in the following scenario. We consider an 
uniformly distributed heavy quark cloud in this hot plasma. We then ask: how does 
the drag force on an external quark change with the density of the quark cloud? 
Like other drag force calculations, we compute it holographically. First, we construct 
the gravity dual of the quark cloud. In the bulk, this represents a black hole in the 
presence of a string cloud. These strings, assumed to be non-interacting, are aligned 
along the radial direction of the bulk geometry and are distributed uniformly over 
R^. We then study the thermodynamics of black hole and verify the thermodynamical 
stability. We also analyze it's gravitational stability under vector and tensor pertur- 
bation in a gauge invariant way. Subsequently, we introduce a probe string with an 
end on the boundary and calculate the dissipative force on the heavy quark following 
the usual approach. At this point, though, the relevance of this work in the light 
of recent quarck-gluon-plasma experiment is not immediately obvious. Nonetheless, 
we would like to make the following comment. The dynamics of a heavy quark (say 
charm) passing through the plasma is usually described by considering it's interaction 
with the medium and the resulting energy loss is calculated. In such calculations, any 
possible effects of other heavy quark due to the back-reaction of the plasma are ne- 
glected. In the context of A/" = 4 SYM, our work can perhaps serve as an attempt to 
compute such back-reaction effects. Within the gauge/gravity correspondence, such 
effects can be modeled in terms of the deformation of the geometry due to finite 
density string cloud. This work shows that the back-reacted gravity background is 
explicitly computable. There are several works, starting with [7], where drag force 
on an external quark has been calculated by introducing D7-brane as a probe in D3 
background. External quark on the D7-brane comes from the end point of an open 
string stretching between the D7 and the horizon of the D3-brane. The usual probe 
approximation here is justified because the free energy of the external quark goes 
as 0{Nc) whereas the plasma, being in adjoint representation, contributes 0{N^) to 
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the free energy. So, in this sense, in the large colour limit, with Nc — t- oo, external 
quark can be treated as a probe. However, when large number of external quarks 
are introduced, the background geometry may get modified as our work indicates. 
Before we go into our computations, here is a note of caution. We study the motion 
of an external quark in an uniformly dense external quark cloud. However, in general, 
this motion is expected to back-react and change the uniformity. Such effects will be 
neglected in this work and, therefore, it is within this approximation that our results 
should be interpreted. 

Gravity dual for external quark cloud 

We consider the (n + 1) dimensional gravitational action given by 

^ - „n+l 



_ „ , dx"^'V^iR-2A) + Sm, (1) 
where Sm represents the matter part of the action. We represent the matter part as 
Sm = -\Y^%j d'^iV^h^^d^X^^dpX'^g^,, (2) 

where we considered g^'^ and /i"^ are the space-time metric and world-sheet metric 
respectively with /i, v represents space-time directions and a, /3 stands for world sheet 
coordinates. 5*^ is a sum over all the string contributions with z'th string having a 
tension 71. The integration in ([2]) is over the two dimensional string coordinates. 
Varying this action with respect to the space-time metric leads to 

R^lv — I^RQiiu + J^g^iu = 87lGn+lTf,u, (3) 

with 

T/^- = _ V 7- / d'^^^J\h^\h-''d^X>^dpX''6'l+\x - X). (4) 

In the above, the delta function represents the source divergences due to the presence 
of the strings. In the following, we will consider the space-time metric of the form 

ds^ = gtt{r)dt^ + grr{r)dr'^ + r"^ 5 abdx"" dx^ ^ (5) 

where (a, h) run over n — 1 space directions. We will further consider strings with 
uniform tensions T and use the static gauge t = C,^,r = The non vanishing 
components of T^'^, following from (jl]), are 

T"--^, T---^. (6) 
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Here we have assumed that the strings are uniformly distributed over n — 1 directions 
such that the density i^ 



aix] 



T J2 - Xi), with a > 0. (7) 



For negative a, T^^ will cease to satisfy the weak and the dominant energy condition^. 
We look for a solution of in AdS space and parametrize the metric acccordingly 
treating a as a constanl0, 

ds^ = -V(r)dt'^ + — — + r'^hiidx'dx^. (9) 
V{r) 

Here hij is the metric on the — dimensional boundary. As for the matter part 
we will focus on to the string cloud for which the nonzero T^i, components ar^ given 
by 



a 



TI=T: = -^^, with a>0. (10) 



r 



The solution which satisfy the Einstein's equation can be easily constructed. It is 
given bj{§ 

nr) = '<+p-—,- („_i),„-3 - (11) 

Here K = 0, 1, —1 depending on whether the (n — 1) dimensional boundary is flat, 
spherical or hyperbolic respectively, having curvature {n — l)(n — 2)K and volume 
Vn-i- In writing down V{r) we have also parametrized cosmological constant as 
A = —n{n — l)/(2/^). With equation ( ITTj) . the metric (|9]) represents a black hole 
with singularity at r = and the horizon is located at V{r) = 0. The horizon has 
a topology of fiat, spherical or hyperbolic depending on the value of K. However, 
our interest in this work, lies in the K = case. In this case of fiat horizon, the 
integration constant m is related to the ADM (M) mass of the black hole as follows, 

M = ^ " ' • (12) 



^To define this properly, we need to think of an IR cutoff in n — 1 directions. 
^For earlier discussions on string cloud/fluid models see [TO l fTT | [T2 | . 

^ Clearly a in ([7]) depends on x. However, in equations PII|) and in ([TT|). a is treated as 
constant. To do this, we have replaced a{x) by an average density as 

a = -^/a(x)d"-^x=-^f]/5r^)(x-X0d"-ia: = -^£ (8) 

Vn-l J Vn-l , J Vn-l , Vn-1 

2 — 1 l — l 

Here, Vn-i is the volume in n — 1 dimensional space after imposing an IR cut-off. Now we consider 
the limit Vn-i going to infinity along with the number of strings N, keeping N/Vn-i constant. 

^It turns out that replacement of Sij by hij in ([5]) keep the components of the stress-tensor same. 

®This is a slight generalization of the metric in |13) . 
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The horizon radius, denoted by r+, satisfies the following equation 



ri 2m 2a 



This allow us to write m in terms of horizon radius as 



0. (13) 



(n - l)r" - 2a/V+ 
The temperature of the black hole is given by 

Note that the zero mass black hole has a non-zero temperature and is given by 

r„ = ^f!izI)Sf. (16) 



2tt V 2aP 

The black hole temperature increases with the horizon size and for large r+, it behaves 
as T ~ r+//^. The entropy is defined as 

S = J T-^dM, (17) 

leading to the entropy density 

Note that s is finite even for black hole with zero mass. The specific heat associated 
with the black hole is 

^ Vn^ijn - l)rl-\n{n - l)rl - 2aPr+) 
dT 4Gn+i{n{n-l)rl + 2{n-2)aPr+) ' ^ ^ 

Now we have a detail look at the thermodynamic quantities just evaluated. First 
of all, if we restrict the temperature to be non-negative, the black hole can have 
minimum radius ^ 

2(ll \ n-l 



It can easily be checked that if we focus on to the region T > 0, there is a single 
positive real solution of ( !T3|) . We also note from ( !20!) and f|T4|) that the mass becomes 
negative at zero temperature 

n 



'^Due to the nature of a dependent term in (jlip. our definition of ADM mass is perhaps ambiguous. 
However, with this definition, entropy of the black hole comes out to be one quarter of the horizon 
area, provided we assume that the equation dS = MdT holds for this black hole. 
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This is somewhat similar to the AdS-Schwarzschild with negative curvature horizon 
[Ti] . We note that, for mass m > m™™, the specific heat f|T9|) continues to be positive 
and is continuous as a function of r^. This suggests thermodynamical stabihty of the 
black hole. Finally, we write down the free energy of this black hole 

(n-l)rl + 2af(n-2)r, , , 
16.(n-l) ■ 

Before we go on to analyze gravitational stability of the black hole, we would like 
to make the following comment. Quite naturally, one may wonder if this black hole 
has a higher dimensional origin. In particular, can this arise, in some near horizon 
limit, from some brane configuration in ten or eleven dimensions after compactifying 
on spheres (with the cloud smeared over the compact manifold)? We indeed tried to 
get this from some bound state configurations of D-branes and strings but have not 
succeeded yet. 

Stability of the fiat black hole 

We now study the stability of the K = ^ black hole geometry using the grav- 
itational perturbation in a gauge invariant way [161 Il3 [El l20l [21]. We consider 
perturbation on a background space time M^"*"^ 

M^+P = A/'2 X KF, (23) 

where the space time metric is, 

ds^ = -f{r)dt^ + + r^Sijdx'^dx^, 

jy^>~~j2~ j.p-1 ~ pj.p-2 ■ V / 

We identify N"^ as a two dimensional space time coordinatized by t and r, whereas /C^ 
is a p dimensional maximally symmetric space coordinatized by Each perturbed 
tensor realized on /C^ can be grouped into scalar, vector, and tensor components such 
that Einstein equations of motion respect the decomposition. Here we do stability 
analysis for tensor and vector perturbations. Scalar perturbation is somewhat more 
involved and will be reported else where. 

Tensor perturbation 

In the case of the tensor perturbation, the metric tensor and energy momentum 
tensor become decomposed in scalar, vector, tensor part with respect to /C^ in the 
following manner [T7] . 

hab = 0,hai = 0, hij = 2r^HTTij 
6T,, = 0,6T^^ = 0,6T; = ttT), (25) 



5 



Tij is the tensor harmonic function defined on JC^. It satisfies the following properties, 

(A + k^)Tij = 
= 0, DjT- = 0. (26) 

Here we note that in /C^ space, A and Dj are realized as the Laplace-Beltrami self- 
adjoint operator and the covariant derivative respectively. For K = 0, can take 
non-negative real continuous values [19]. Gauge invariant quantities like Ht and tt 
are function of variables belong to M"^ spacetime. [2D]- 

Now substituting all the variations in the perturbed Einstein equation, we get the 
master equation of tensor perturbation [17]. 

T 2 

DHt + -Dr- DHt - -^Ht = -k^tt (27) 
We introduce a new variable 

$ = t'p/'^Ht, (28) 
and substitute it into the master equation. It takes following canonical form, 

- ^ = 0, (29) 



where Vt is defined as. 



V,= Llkr' + ?^ + ?^^]. (30) 



According to (jl]) the energy momentum tensor is constructed with the spacetime 
vector X^{t, r) which does not contribute to the linear order of gauge invariant tensor 
perturbation. Therefore we set tt to be zero in fl2I?l) . It is clear in plot-1 of figure ^ 
that for higher dimensions Vt is always positive beyond horizon. So the black hole 
geometry is stable under tensor perturbation. 
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Vector perturbation 

In the case of vector perturbation, metric and energy momentum tensor are de- 
composed in terms of vector harmonics Vi as well as vector harmonic tensor l^j|17]. 

hab = 0, hai = rfaVi, hij = 2r^HTVij 
6Tat = 0, 5Tt = rr'^V,, 5T; = vrVj (31) 

The vector harmonics are defined as 

(A + kv^)Vi = 0, DiV' = (32) 

/^From vector harmonics we can construct vector type harmonic tensor, 

V,, = -^iD,V,+D,V^, 

(A + 4)\/,, = 0, Vi = 0, D,V^ = ^V,. (33) 
The gauge invariant parameters for K = are given by 

Fa = fa + rDa{^),TT^T\ (34) 
Ky 

Upon considering the perturbations of the Einstein equation and the conservation 
law of energy momentum tensor, master equation arising from the gravitational per- 
turbation with the source term takes the following form [19], 



r^D,{-D^n) - = -^rV^D„(rr,) (35) 



where. 



^P-ipa ^ ^abjj^^ ^ rP-^^r\ (36) 

ky 

Now introducing the change of variable 

$ = r-P/^Q (37) 

we recast the master equation in a canonical form, where the effective potential for 
vector perturbation comes out as [21], 



V.,= /|,/ + I<P±^;_|;T (38) 

The plot-2 in figure (l2|) implicates that beyond horizon, Vy is not always non- 
negative for p > 3. We follow S — deformation method [T7] to construct modified 
effective potential 

VV = + - S\ (39) 




(a) Plot-1 



(b) Plot-2 



Figure 1: Plot 1 shows, for various dimensions, the effective potential Vr in tensor 
perturbation is positive beyond horizon radius . Plot 2 shows the effective potential Vy 
in vector perturbation is not always non-negative for p > 3. In both cases horizontal 
axis is normalized with respect to black hole horizon radius r+. 

where S is an arbitrary function of r. If we choose S = we get the modified 
effective potential. 



furnishes the sufficient condition for the stability of the black hole. 

Having constructed this black hole geometry we compute the drag force on an external 

quark moving in external quark cloud 

Dissipative force on an external quark moving in the heavy quark cloud 

We now like to calculate the dissipative force experienced by the external heavy 
quark moving in the cloud of heavy quarks. Our aim is to study the force as a function 
of the cloud density. Calculational procedure to evaluate drag force on an external 
quark is by now standard. This can be found, for example, in [3 [HI |9] . We will follow 
the notations of [9j. The drag force on a very massive quark with fundamental SU{N) 
charge at finite temperature is calculated holographically by studying the motion of 
a string whose end point is on the boundary. This end point represents the massive 
quark whose mass is proportional to the length of the string. We will consider here 
the gauge theory on coordinatized by x^,x'^,x^. This means, for the purpose of 
this computation, we only consider K = 0,n = 4 case of the black holes discussed 
previously. 

Let us consider the motion of a string only in one direction, say x^. In static 
gauge, t = ^°,r = the embedding of the world-sheet is given by the function 
r). The induced action of the string in our case follows from a straightforward 
computation 




S 



1 



I 



dtdr A / 1 + 



3r^ — 2aPr — GmP 



{drx^y 



3r4 



{dtx^y 



2na' 



3/4 



3r4 — 2aPr — 6mP 



(41) 



8 



where we have scaled by /. 

The ansatz that describes the behaviour of the string with attached quark moving 
with constant speed v along is given by [9] 

x^{r,t) = vt + i{r), (42) 

for which fj^Tl) simplifies to 



S = / dtdr\ 1 + (drO^ -v^. 43 

The momentum conjugate to ^(r) is 

1 (3r4 - 2afr - 6mP)dr^ , ^ ^ , 

TTf = — =. (44) 

Equation of motion can be obtained by inverting this equation for dr^. However, as 
in ^Qj, to get a real ^, the constant of motion tt^ has to be set to 

where r„ is the real positive solution of the equation 

3(1 - v'^yl - 2a/V„ - 6m/2 = 0. (46) 

Though this equation can be solved explicitly, the solutions are not very illuminating. 
However, it is easy to check that there is only one real positive solution. Substituting 
this solution of r„ in we get tt^. The dissipative force is then given by [0] 

1 (47) 



27ra' P ' 



with r„ given by the positive real solution of f l46|) . Now we wish to rewrite the 
expression of the dissipative force in terms of gauge theory parameters. Along this 
line, we solve ( fT5|) for r+, 

r+ = ^A(T,6), (48) 
where h is the scaled quark cloud density, h = a/l^ and A is given by, 

A{T,b) = [{2(96 + Att^T^) + 6^/6(96 + Stt^T^)}^^^ + 

92/3 2^ 

2nT{l + = }] (49) 
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-50^ • |- 

(a) Plot 1 (b) Plot 2 

Figure 2: Plot 1 shows the variation of F as a function of T for the values of quark 
density 6 = (solid), 0.5 (dashed), 2 (dotted) respectively. Plot 2 shows the variation 
of F as a function of b for the values of T = .1 (solid), .5 (dashed), 1 (dotted) 
respectively. We see in both cases the larger the quark cloud density as well as 
temperature, the more is the dissipative force. 



Substituting (H8l) and the following useful relation 



— = gYM^N, (50) 
a' 



in the expression of the dissipative force fl47j) . we get the modified form 

"^'■^/^^v^. (51) 



727r ^ r. 



Here gyM is the Yang-Mills(YM) gauge coupling and is the order of the gauge 
group SU{N). We are able to solve the ratio r^'^/r^'^ in a closed form by substituting 
( |T3l) into f H6l) . The relevant equation takes the following form, 

(1 _ v^){^) lii^((!:!i) - 1) _ 1 = 0, (52) 

^ (A(T,6))'^V+^ ' ^ ^ 

It turns out that the real positive solution of fl52|) is expressible in terms of A{T, b) 
and b itself. Denoting the solution as f{A, b) and plugging it back into flSTj) we achieve 
the form of dissipative force expressible in terms of gauge theory parameters 

A' 



727r 



^/^^^vfiA,by. (53) 



We note here that f{A, b) is an explicitly computable function. 

We would now study ( 153|) for different values of heavy quark density and for fixed 
T. As for an example, it is interesting to check that if the temperature is fixed 
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at the value Tq as mentioned in f|T6|) the dissipative force behaves as F ~ 

_52/3 

where b is now the density of the quark cloud. Also for T = 0, A{T, b) in poj) 
simplifies significantly resulting the dissipative force to vary as F ~ —6^/^. For 
generic temperature and small b, it is possible to have a power series solution of flS^ 
in 6. However, for appreciable density, we find it more suitable to analyze F in terms 
of plots. In figure ([2]) plot 1 shows the behaviour of the drag force as a function of T 
for different b. For fixed T, we clearly see that the force becomes stronger with the 
quark densitjl^ 

Conclusion 

In this paper, using the AdS/CFT duality, we have computed the dissipative force 
experienced by an external heavy quark with fundamental SU{N) charge moving in 
the heavy quark cloud at finite temperature. In the dual theory, we have considered 
motion of a string (with one of it's end point at boundary) in a + 1 dimensional 
background with flat boundary pervaded with string cloud as the matter source. 

The geometry implies existence of a black hole parametrized by it's mass m and 
the string cloud density a. The black hole turns out to be thermodynamically stable 
and it resembles AdS-Schwarzschild black hole with negative curvature horizon. We 
have been able to check the gravitational stability of the geometry for tensor and 
vector perturbations. 

In the above scenario, we have computed the drag force exerted on the external 
quark. The most general form of the dissipative force turns out to be a complicated 
function of temperature of the boundary theory T and the re-scaled quark cloud 
density b. However for the temperature corresponding to massless black hole in the 

2 

dual gravity theory, it behaves like F ~ —63. We have plotted it with respect to both 
T and b separately while keeping one of them constant at a time. Both plots exhibit 
an enhancement in the drag force in the presence of evenly distributed quark cloud. 

There are few issues that we hope to look into. Firstly, what is the higher di- 
mensional brane geometry whose near horizon geometry contains the black hole that 
we are considering? Secondly, is the gravitational background stable under scalar 
perturbation? 
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